Log-periodic oscillations for diffusion on self-similar finitely ramified structures 



L. PadillaQ H. 0. MartinQ and J. L. Iguainf] 

Instituto de Investigaciones Fisicas de Mar del Plata (IFIMAR) and Departamento de Fisica FCEyN, 
Universidad Nacional de Mar del Plata, Dean Funes 3350, 7600 Mar del Plata, Argentina 

Under certain circumstances, the time behavior of a random walk is modulated by logarithmic 
periodic oscillations. The goal of this paper is to present a simple and pedagogical explanation of the 
origin of this modulation for diffusion on a substrate with two properties: self-similarity and finite 
ramification order. On these media, the time dependence of the mean-square displacement shows 
log-periodic modulations around a leading power law, which can be understood on the base of a 
hierarchical set of diffusion constants. Both the random walk exponent and the period of oscillations 
are analytically obtained for a pair of examples, one fractal, the other non-fractal, and confirmed 
by Monte Carlo simulations. 
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I. INTRODUCTION 

Diffusion on non-Euclidean media is not a new re- 
search area. For example, the underlying mechanisms 
of the sub-diffusive behavior, characteristic of fractals, 
were discovered some decades ago. Indeed, it is by now 
well established that, on these objects, the spreading of 
the probability density function is continuously retarded 
because of the presence of holes of all sizes (the inter- 
ested reader can refer to [IH5], and references therein). 
However, in the last years, considerable effort has been 
dedicated to investigate an outstanding phenomenon. It 
has been repeatedly reported that, on some determinis- 
tic fractals or self-similar graphs, the time behavior of a 
random walk (RW) is modulated by logarithmic-periodic 
oscillations [BHTTj. Similar effects were observed for bi- 
ased diffusion of tracers on random systems |12H15j , and 
out of the domain of diffusive motion, examples have 
been detected in earthquakes [THl [H], escape probabil- 
ities in chaotic maps close to crisis |18) . kinetic and dy- 
namic processes on random quenched and fractal media 
[TM2"2"] . diffusion-limited aggregates [33], growth models 
[24], and stock markets near a financial crash [25ll28j . 
There is general agreement that these oscillations are a 
manifestation of an inherent self-similarity [5^, respon- 
sible for a discrete scale invariance [50] . 

A wide class of systems exhibiting log-periodicity is 
that of self-similar finitely ramified structures. Indeed, 
several researchers have found that, on these substrates, 
the root-mean-square displacement (RMSD) of a single 
RW as a function of time, follows an anomalous power 
law modulated by logarithmic periodic oscillations [9HTT] 
(for related behavior, see [BHH] and [IS])- In this paper, 
we revisit this problem, with the purpose of explaining 
the origin of the behavior described above in a simple 



and comprehensible way. The work is a natural continu- 
ation of a previous one [31] , where we studied log-periodic 
modulation in one-dimensional RW. 

Our approach is two-fold: theoretical analysis and nu- 
merical simulations. In the theoretical part, we treat the 
diffusion problem in a quite general way. The analysis is 
facilitated by the fact that a part of a finitely ramified 
object can be separated from the rest by cutting a finite 
number of connections, which, because of self-similarity, 
does not depend on the size of the part to remove. For 
pedagogical reasons, the procedure is illustrated by two 
simple examples but the method should be easily gen- 
eralized for other finitely ramified substrates. The first 
model of substrate, denoted I, is the well-known Vicsek's 
snowflake fractal structure [35]. The second, Model II, is 
also finitely ramified but may be considered a trivial frac- 
tal (with a fractal dimension of 2 and without holes of all 
sizes) . Model I is an example of what we will call a perfect 
diffusive self-similar structure; Model II corresponds to 
what we will call an asymptotically diffusive self-similar 
one. For both models we predict a time-modulated os- 
cillatory behavior, and calculate the RW exponent and 
the period of the oscillations. The goal of the numerical 
part is to get some independent confirmation of the de- 
rived properties. We study, using standard Monte Carlo 
(MC) simulations, the time evolution of a particle diffus- 
ing on each substrate, and compare the findings with our 
theoretical predictions. 

The paper is organized as follows. In Sec. [IT] we define 
the models, discuss the causes for modulated power-law 
behavior and derive the RW time evolution. The outcome 
of MC simulations are presented in Sec. |III| Finally, we 



give our conclusions in Sec. IV 



II. ANALYTIC APPROACH 
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The above mentioned substrates are built in stages on 
a square lattice and the result of every stage is called a 
generation. The structure of each generation is a peri- 
odic array of basic or unit cells, which consist of sites 
connected by bonds. We denote by L the linear size of 
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FIG. 1: The unit cells of Model I. The zeroth, first, and second 
generations are shown in (a), (b), and (c), respectively. L (3 
in this case) is the linear size of the first generation. 



the unit cell that corresponds to the first generation. 

On any generation, the motion of a single particle, ini- 
tially located at a symmetry point of the structure, occurs 
stochastically. The particle jumps, with a hopping rate 
k, only between NN substrate sites which are connected 
by a bond. 



A. Model I 



(A 2 x)(t) = 2D&>t, (1) 

which holds for a time t longer that the average time for 
the particle to escape from the initial unit cell. Because 
of the x o y symmetry, the same time dependence holds 
for the mean-square displacement in y direction. 

It is not hard to convince oneself that the zeroth- 
generation substrate is the simple square lattice (see 
Fig. [l]-(a)), and then 

D<® = k. (2) 

The first-generation substrate is shown at the top of 
Fig. [2j However, regarding x-direction diffusion, the 
whole substrate and the string of cells displayed at the 
bottom of Fig. [2] with periodic boundary conditions in 
the y-direction, lead to equivalent problems. We exploit 
this equivalence and calculate the diffusion coefficient of 
that one-dimensional array of equivalent cells. Following 
the steady-state method of Ref . [33] , we obtain 

D« = (3/5)fc. (3) 

The reduction to a quasi one-dimensional problem is 
possible for every n, and, after some simple algebra, it is 
found that 



The building process is illustrated in Fig. [I] Parts (a), 
(b) and (c) of this figure show the unit cell for the zeroth, 
first and second generation, respectively. It is easy to see 
that, for this model, L = 3a, where a is the distance 
between nearest-neighbor (NN) sites (a = 1 in the rest of 
the article). It is also apparent from this figure that the 
second-generation basic cell has a linear size L 2 and is 
built from the first-generation one in a self-similar way. 
In general, the n-th generation basic cell has a linear size 
L n and, for n > 1, a central part of a linear size 

L (n-rn) 

(with m — 1,2, ...,n — 1) can be separated from the rest 
by cutting only four bonds. 

For the n-th generation, a two-dimensional periodic 
substrate is built by connecting the corresponding basic 
cell (the first-generation substrate can be observed at the 
top of Fig. [2]). The final full self-similar substrate, we 
are interested in, is obtained after an infinite number 
of iterations. Note that, in this limit, the unit cell is 
identical to the snowflake fractal. 

We proceed now to analyze the behavior of the dif- 
fusing particle. It is useful to remember that, on any 
periodic substrate, normal diffusion should be observed 
if time is long enough for the RW to be influenced by the 
structure periodicity. Thus, for the n-th generation sub- 
strate, a diffusion coefficient can be defined through 
the time dependence of the mean-square displacement 
(A 2 x)(t) = ([x(t) - x(0)} 2 ) in the x direction, 



D {n) = (3/5)™fc, for n = 0,1, 2,.... (4) 

On the base of this result, we can anticipate that, 
on the full self-similar structure, a RW will show a log- 
periodic modulated behavior. The key point here is that 
the diffusion coefficients satisfy 

D(")/B ( " +1) =6, for n = 0,l, 2,..., (5) 
with S = 5/3. 

Let us start noting that, since at a time t the par- 
ticle will be located in a region of linear size of the 
order of the RMSD */ (A 2 x) (t), for a time such that 
AL n < y / (A 2 x)(t) < AL n+1 (A is a constant of the order 
of one) , it will be impossible for the RW to distinguish the 
full self-similar substrate from the (n + l)-th-generation 
one. Everything will thus happen as on the latter and 
one should expect that 

(A 2 x)(t) = 2D (n) t, (6) 

for t in that interval. 

At later times, when the RMSD is of the order of L n+1 , 
the particle will start to diffuse as if located on the (n+1)- 
th generation substrate. The relation 
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FIG. 2: First generation of Model I. Top: the substrate built 
with the basic cell shown in Fig. [T]-(b). Bottom: the infinite 
one-dimensional string of cells used to compute the diffusion 
coefficient D w (see Fig.2 in Ref. [53]). The arrows indicate 
periodic boundary condition in the y direction meaning that 
when, for example, the RW at site a (b) jumps upward (down- 
ward), it will arrive at site b (a). 



(A 2 x)(t) = 2D ( - n+1 h 



(7) 



will thus hold for AL n+1 < (A 2 x)(t) < AL n+2 . 

Because of the substrate full self-similarity, we should 
observe a sequence of length-dependent diffusion coeffi- 
cients, giving rise to the qualitative behavior depicted in 
Fig. [J This is a sketch of the mean-square displacement 
as a function of time (thick curve), which has a power 
law form modulated by a log-periodic amplitude. That 



(A 2 x)(t) = Ct 2 »f(t), 



(8) 



where v is the RW exponent, f(t) is a log-periodic func- 
tion which satisfies /(tr) = f(t), with the logarithmic 
period log(r), and the constant C is obtained by asking 
that the log-time average of log(/) over one period be 
zero (see, for more details, Fig. [5]). 



FIG. 3: (Color online) Schematic of the mean-square dis- 
placement as a function of the time, shown by the thick 
red curve. The length of the segment be is log(2D^) — 
log(2D (n+1) ) = log(<5), because of Eq. J5I. From the slopes 
(= I) of the full straight lines (representing the normal dif- 
fusion behaviors, {A 2 x} = 2D (n) t), one gets that the seg- 
ments ad and cd have the same length or, equivalently, that 
log(r) = log(L 2 ) + log(5). The dashed straight line represents 
the global power law (A 2 x) ~ i 2 ", with 2v = log(i 2 )/ log(r). 
More details in the text. 



In Fig. [3j two groups of inclined straight lines were 
drawn as a guide. The dashed line has a slope of 2v and 
corresponds to the (A 2 :r) global power-law trend. The 
solid lines have slopes of 1 and represent normal diffusion 
in each of the different generation structures. Both r and 
v can be expressed in terms of L and 5. It is clear from the 
slopes of the straight lines that log(r) = log(L 2 ) + log(<5) 
(solid line) and 2v = log(L 2 )/ log(r) (dashed line), which 
is equivalent to 



5L 2 = L 1/V 



and 



logi 



(9) 



(10) 



Note that, since 5 > 1, anomalous sub-diffusion results 
(y < 1/2, see EqfTo]), and that, according to the sketch 
in Fig. [3j the amplitude of the modulation increases with 
the increase of S. 

We would like to stress again that the self-similarity 
in the mean-square displacement, schematically shown 
in Fig. [3] and mathematically described by Eq. pi), is a 
direct consequence of It is because of the latter 

that, in Fig. [3j the distance between any pair of near- 
est solid straight lines is a constant, and that two nearest 
equivalent points (like a and b) are always related by 
the transformation (t -> rt , (A 2 x) -> t 2v (A 2 x)). The 
set of equations §5§ also allows us to obtain analytically 
both the random walk exponent and the period of the 
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FIG. 4: The unit cells of Model II. The zeroth, first, and sec- 
ond generations are shown in (a), (b), and (c), respectively. 
L — 3, as in Model I. The dashed lines correspond to the 
bonds which are removed from the immediately previous gen- 
eration. Note that, for every order, two NN unit cells are 
connected by only one bond. 



oscillations. We will refer to Model I as a perfect diffu- 
sive self-similar structure. Other examples of this kind 
of structures are the fractals shown in Fig.(l) of Ref. [PT] . 



B. 



Model II 



For this model, the basic cells for the zeroth, first and 
second generation are shown in Fig. [4] The full self- 
similar substrate is obtained when the generation order 
n goes to infinity. As for Model I, the width of the n-th- 
generation unit cell is L n (where L = 3 is the width of 
the first-generation one), and, for n > 1, a central part 
of linear size L n ~ m (for m < n) can be separated from 
the rest by cutting only four bonds. Again, this allows 
us to reformulate the RW problem on a one-dimensional 
array (similar to that shown in Fig. |2]-bottom), and to 
compute the diffusion coefficients following the steady- 
state method [33] . 

It is evident that L> (0) is given by Eq. |2j) . For the next 
two generations we obtain 



and 



DW = (l/2)fc, 



= (45/174)fc. 



(11) 



(12) 



We stop here because the calculations become more te- 
dious for higher orders. However, at this point, we are 
already able to grasp the general trend, and to make a 
comparison between the two models. 

From the values of the first few diffusion coefficients 
(Eqs. ([2J), pdj ) and (12)), we immediately observe that 
the ratio 1W>/D^ = 2 is different than D^/D^ = 
87/45 ~ 1.933. In general, an analysis of the structures 
indicates that, for Model II, the equations ^ should be 



D {n) /D {n 



+i) 



(13) 



where, instead of a constant 8, we have now a sequence 
8 n . In spite of this modification, the qualitative behavior 
shown in Fig. [3] still holds, but the exponent v should be 
replaced by a sequence 



1 

2 + ^ 



and the constant period r by 

< n — *-* 



(14) 



(15) 



It is interesting to note that v n and r„ can naturally be 
interpreted as a local RW exponent and a local period, 
respectively. 

Substituting the parameters above (<5o = 2 and <5i = 
87/45) into Eqs. (fit]) and (fij>]), one obtains v ~ 0.3801, 
vx ~ 0.3846, log(7u) ~ 1.2553 and log(n) ~ 1.2406. 
Even if we do not calculate 8 n for every n, in the next 
paragraph we argue that an oscillating modulation is al- 
ways present, and that a log-periodic behavior appears 
in the long-time limit. 

Since the above results seem to indicate that 8 n is a de- 
creasing sequence, a crucial question is whether, for large 
enough n, it approaches the trivial limit 8^ = 1, equiv- 
alent to normal diffusion (see (Eq. [m])). To see that this 
is not the case, look at Fig. [3] and observe that whenever 
the distance in the x direction from the center of symme- 
try exceeds a threshold L n /2, the number of connections 
per unit length is suddenly reduced by a factor of 3. It 
will be increasingly difficult for a random walker to reach 
large distances, and a decrease in the diffusion coefficient 
will happen at each of these characteristic lengths; thus 
8 n > 1 should hold for every n. Moreover, a compar- 
ative analysis of the structures composing Model II al- 
lows us to obtain a better (though rather crude) lower 
bound. Indeed, it is possible to show that the inequality 
8 n > 7/5 is always true. We can thus conclude that, as 
n goes to infinity, the sequences 8 n , v n and t„ have well- 
defined limits (5oo > 7/5, i^oo < 0.4336 and r m > 12.6, 
respectively), and that an asymptotic regime emerge, in 
which an overall sub-diffusive behavior is modulated by 
log-periodic oscillations. We then say that Model II is an 
asymptotically diffusive self-similar structure. 

It is worth to mention that in this full self-similar struc- 
ture the distribution of sites is that of the square lattice. 
Therefore, in the sense of that distribution, this substrate 
is a trivial fractal of dimension d m = 2 and without holes 
of all sizes, which could be considered the cause of sub- 
diffusion. One may ask if, regarding RW, the distribution 
of bonds is more important than that of sites. However, 
the related bond fractal dimension db is also db = 2: as 
the counting-box linear size is increased by a factor of L 
(from L n to L n+1 ), the number of bonds increases from 
N b to L 2 N b — 4 (see Fig. [I| for n = 1), which is equiv- 
alent to a L 2 increase factor, for large enough n. Thus, 
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Model II exhibits logarithmic time periodicity without a 
(nontrivial) fractal structure. 



III. 



NUMERICAL RESULTS 



As mentioned in Sec. |TJ in this section, we explore dif- 
fusion using MC simulations according to models I and 
II. For each model, we perform standard MC simulations 
of a single RW on the n-th generation basic cell. Every 
RW starts at the center of symmetry of the cell. The 
value of n is always chosen large enough to prevent the 
RWs from reaching the cell borders during the simula- 
tion. Working on this cell is thus equivalent to working 
with the full self-similar structure. In all simulations the 
time step At is set to 1 and the hopping rate k to 1/4. 
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FIG. 6: Log-log plot of the mean-square displacement versus 
time for Model II. The straight line has a numerically fitted 
slope 2v = 0.77. The inset is a plot log( A 2 x}/Ct 2 " for the 
same data and the curvilinear line was obtained as in Fig. [5] 
The period r is here obtained from Eq. [9] with the numerical 
value of v {vfu)- 



From these figures it is indeed clear that the predictions 
in Eqs. (pi), (10 1, (14 1 and (15) are consistent with the 
numerical findings. Note also that the amplitude of the 
oscillations is larger in the case of Model II, consistent 
with its higher value of S. 



IV. 



CONCLUSIONS AND DISCUSSION 



FIG. 5: The mean-square displacement (A 2 x) as a function 
of time for Model I. The straight line has a slope 2v, with 
v ~ 0.4057 obtained from Eq. (10 1. The inset shows a plot 



of log(A x)/Ct " vs. logt for the same data, where C is a 
properly chosen constant. The curvilinear line represents the 
first-harmonic approximation A sin(27r log t) / log(r) +«). The 
period r is given by Eq.(|9|. A and a are fitted constants. 

The numerical results for the mean-square displace- 
ment are plotted as a function of time in Figs. [5] and [6] for 
Model I and Model II, respectively. It is apparent from 
these figures that (A 2 x)(t) satisfies a modulated power 
law. The modulations are more easily observed in each 
inset, where \og({ A 2 x) / Ct 2v ) versus log(f), is plotted us- 
ing the same data as in the corresponding main plot (C is 
a constant chosen to have the oscillations centered around 
zero). In the case of Model I, the value of the RW expo- 
nent (y = 0.4057) was computed from Eq. (10), while for 
Model II it was numerically fitted (vfu = 0.385). Note 
that, in the last case, the value is very close to the an- 
alytical approximation v\ — 0.3846, indicating that the 
convergence of the sequence is quite fast. The curvilin- 
ear lines are of the form Asin(27r log(i)/ log(r) +a), i. e., 
the first-harmonic approximation of a periodic function 
with period log(r), where A and a are fitted parameters. 



We have studied the time evolution of a single RW on 
a self-similar finitely ramified structure. The arguments 
employed in this work show that the time-modulated os- 
cillatory behavior these systems exhibit originates from 
finite ramification. 

Schematically, at a length L n , the characteristic unit 
cell may be visualized ge with a given number of 

gates. While, inside the cage, the RW evolves normally, 
the situation is significantly different in the border, where 
the diffusive motion is effectively slowed down due to the 
small number of gates. A relevant parameter is the ratio 
between the number of gates and the cage perimeter. Be- 
cause of finite ramification and self-similarity, this ratio 
decreases as n increases, and thus, the larger the cage, 
the harder it will be for the RW to find a way to get out. 

The RW dynamics can thus be interpreted as a se- 
quence of normal diffusion behaviors, each characterized 
by a diffusion coefficient, which decreases as the length- 
scale increases. The oscillatory modulation amplitude 
of the mean-square displacement is thus an emergent 
property of the transitions between two normal diffusion 
regimes. Because of the substrate self-similarity, the os- 
cillatory behavior survive in the long-time limit. 

We have shown that there exist two wide classes of 
self-similar finitely ramified structures. For pedagogical 
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reasons, we have illustrated their properties by studying 
two models: the snownake fractal (Model I) and a (triv- 
ial) fractal with a fractal dimension equal to two (Model 
II). We say that Model I is perfect diffusive self-similar 
because, on this substrate, a single RW follows, at all 
time-scales, a power law modulated by log-periodic os- 
cillations. In contrast, we say that Model II is asymp- 
totically diffusive self-similar since, in this case, the log- 
periodic modulated power law regime is only reached for 
long enough times. For Model I we have obtained an- 
alytically both the global RW exponent and the period 
of the oscillations. Instead, for Model II, we have shown 
that a convergent sequence of local exponents and peri- 
ods exists, and have calculated its first few elements. To 
check the validity of our approach, MC simulations, were 
also carried out for each model. The numerical results 
confirm our theoretical predictions. 

We would like to stress that the (nontrivial) fractal 
character of the substrate is not a necessary condition for 
anomalous sub-diffusion modulated by log-periodic oscil- 



lations. It is a general belief that, in a fractal, a RW be- 
comes sub-diffusive because of the holes of all sizes, which 
hinder access to some regions. However, we have found a 
structure (Model II) without those holes but leading to 
sub-diffusion. 

Let us finally observe that, so far, we have considered 
RWs that always start at the center of symmetry of the 
structure. However, as shown elsewhere [34) . for long 
enough times, the diffusion became independent of the 
initial position, and the qualitative behavior of Fig. [3] 
will hold for L n I, where t is the distance between the 
initial position and that center. 

In summary, for a large set of substrates, we have found 
that the RW sub-diffusive behavior modulated by a log- 
periodic amplitude can be viewed as an emergent prop- 
erty of self-similarity and finite ramification. 

This work was supported by the Universidad Nacional 
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